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Осесимметричный изгиб круглой многослойной пластины на упругом основании 
сложной структуры* 


С. М. Айзикович, С. С. Волков, А. В. Мелконян 


Получено в аналитическом виде приближённое решение задачи об изгибе круглой многослойной пластины 
постоянной толщины, лежащей на упругом основании сложной структуры. Пластина изгибается под действием 
осесимметричной распределённой нагрузки и реакции со стороны основания. Упругое основание представляет 
собой непрерывно-неоднородный по толщине слой (покрытие), лежащий на однородном полупространстве 
(подложке). Модуль Юнга в зоне сопряжения покрытия и подложки имеет существенный скачок. Для пластины 
рассмотрены два случая граничных условий: условия закреплённого и свободного края. Построенное прибли- 
жённое аналитическое решение задачи эффективно в широком диапазоне как геометрических параметров 
(толщина неоднородного слоя и радиус пластины), так и физических параметров (гибкость пластины и упругие 
свойства покрытия и подложки). Методом интегральных преобразований контактная задача сводится к реше- 
нию системы интегро-дифференциальных уравнений. Полученные формулы могут быть использованы для рас- 
чёта характеристик контактного взаимодействия многослойной пластины с основанием сложной структуры в 
зависимости от граничных условий и характера нагрузки на пластину. 
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Введение. Применение функционально-градиентных материалов существенно влияет на все ха- 
рактеристики контактного взаимодействия [1-3]. В случае, когда во взаимодействие вовлечены 
тонкие гибкие элементы (пластины), на перераспределение контактных давлений в зоне контакта 
влияет как неоднородность свойств взаимодействующих элементов, так и жёсткость пластины. Это 
нужно учитывать при расчёте взаимодействия гибких элементов (пластин) с неоднородными струк- 
турами, что приводит к необходимости рассмотрения контактной задачи о взаимодействии пла- 
стины и неоднородного основания. 

Простейшая осесимметричная задача контактного взаимодействия как для слоистого, так и 
непрерывно-неоднородного покрытия упругого полупространства рассмотрена в работе [4]. Задача 
об изгибе пластины на упругом изотропном и однородном основании рассматривалась в работах 
[5, 6]. Решение строилось путём представления контактных напряжений в виде степенного ряда, с 
последующим определением коэффициентов разложения из бесконечной алгебраической системы 
уравнений. 

Методом ортогональных многочленов такая задача решалась в работах [7, 8], а методом 
коллокации по чебышёвским узлам — в работах [9, 10]. При этом возникала необходимость постро- 
ить решение некоторых бесконечных систем линейных алгебраических уравнений и ставилась про- 
блема исследования сходимости полученного решения к точному. В работах [11, 12] для решения 
задачи применялись асимптотические методы типа «больших Л» и специальных ортогональных 
многочленов, что позволило получить основные характеристики решения в нескольких формах, 
каждая из которых эффективна в своей области изменения характерных параметров задачи. 





* Результаты работы получены при выполнении проекта, поддержанного грантом РФФИ № 13-08-90916-мол_ин_нр. 
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Отметим, что большинство известных решений эффективны только для жёстких пластин. И 
очень немногие, в частности, представленные в [11, 12], эффективны или для гибких, или для 
жёстких пластин, каждое в своей области. 

Интерес к решению задачи и её актуальность сохраняется и в настоящее время. Так, в ра- 
боте [13] решение строилось с использованием разложения напряжения в двойной ряд Фурье. Ана- 
логичный подход использовался в работе [14]. Апагеа В. О. 5Й\а с соавторами развил численные 
методы решения задачи [15]. В работе [16] было получено приближённое аналитическое решение 
осесимметричной задачи об изгибе пластины со свободным краем, лежащей на функционально- 
градиентном основании. Для этого использован двусторонне асимптотический метод [17]. 

В настоящей работе метод [17] использован для решения задачи об изгибе многослойной 
пластины. 

Постановка задачи. Круглая многослойная пластина радиуса А и постоянной толщины А лежит 
на поверхности 2 = 0 упругого полупространства, состоящего из неоднородного мягкого слоя (по- 
крытия) толщины Н (-Н <2 <0) и однородного полупространства (подложки), причём /, — по- 
стоянная толщина каждого слоя пластины. С полупространством связана цилиндрическая система 
координат г,ф,2 ; координата готнесена к радиусу пластины К, координата 2 — к толщине покры- 


тия Н, ось 2 нормальна поверхности и проходит через центр пластины. Под действием осесиммет- 
ричной распределённой нагрузки р”(г) и реакции со стороны слоя пластина изгибается. Прогиб 


пластины обозначается функцией и/' (г). 


Коэффициенты Ламе мягкого слоя меняются с глубиной по закону: 
Л.(2), -1<2<0 
кое ео 


^., —©<2<-1, 
| 


М, (2), -15250, 
М,, 0 <2<-1. 
При расчётах для задания характера неоднородности покрытия удобно использовать мо- 
дуль Юнга, а коэффициент Пуассона зафиксировать и считать постоянным. Известно [18], что 
ЗА +2М Л 


Л+М ' = 2(^+М)° 
В качестве параметра, характеризующего отличие упругих свойств слоя от подложки, вве- 
дём величину 


Е =М 


Е, 
Е (1) 


1 
где Е, (2) и Е, — модули Юнга покрытия и подложки соответственно. 


Уравнение изгиба многослойной пластины имеет вид: 
и (г) = р” (г) - 9. (г), 


нии, (г) = 9, (г)-9,, (г), 0<г<1, (2) 


ких (^) = к... (Г)-9" (г), 
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где Ь, =2, 


‚|+ 


в ла 
а? г 


2 
— дифференциальный оператор изгиба АЙ пластины; 4“ (г) — контактные 


напряжения под пластиной; и/, (г) — прогиб каждого слоя; О, — цилиндрическая жёсткость 


Ей пластины. 

Для напряжений, возникающих между слоями пластины, выполнено следующие условие 

сопряжения: 
РЭР)» ФРЕЬ (3) 

Если все слои жёстко сцеплены друг с другом, то и’, (г) =и/, (г) =...=ии‚ (г) =и/" (г). Учи- 
тывая это обстоятельство, а также равенство (3), система уравнений (2) может быть представлена 
в виде: 
а 1а 
ата 

Таким образом, вместо многослойной платины можно рассматривать однослойную пластину 
с цилиндрической жёсткостью, равной сумме цилиндрических жёсткостей каждого слоя в много- 
слойной пластине. 

Если пластина неоднородна по толщине, её модуль Юнга Е’, (2) при этом изменяется по 


рим" (=) -9' (тер = Ур, =| } 0</<1. (4) 


толщине пластины, а коэффициент Пуассона у„ постоянный, то, согласно [19]: 





(2). (5) 


Введём следующую замену переменных: Л = Н/К, и" (г) =и(г’)В, р’(г)= р(г’) Ок, 
4’ (г)=9(г’)Ок?. Знак штриха далее опускаем. 

Рассматривается два случая граничных условий: первый, когда края пластины закреплены, 
в этом случае функция прогиба пластины ихг) удовлетворяет условиям: 


ди д 














=0, Ди =0, 6 
г Г=1 д ) 
и второй случай, когда имеют место условия свободного края: 
2 
ИИ", 0. (АИ 50, (7) 
Ге] г ог га Ге] ты 








где у — коэффициент Пуассона пластины; Л — оператор Лапласа. 
пл 


Методом интегральных преобразований контактная задача сводится к решению системы 
уравнений: 


Е.И (г) =р(г)-9(г), 0<г<1, (8) 
[9 (а) (а^)л (аг) да = зи (г), 0<г<1, 
0 
ь (9) 
[9(а)^ (аг) аа =0, г>1, 
0 
где (и) — трансформанта ядра; $ = ©9АЗО-! — параметр, характеризующий изгибную жёсткость 


пластины. 
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Построение решения. Для однородной пластинки со свободными краем, лежащей на неоднород- 
ном полупространстве, в работе [16] построено решение задачи в аналитическом виде. Учитывая 
(4), (5), это решение можно использовать и для неоднородных пластин. 

Построим решение задачи в аналитическом виде, в случае, когда края пластины закреп- 
лены. Для этого представим функцию прогибов в виде ряда по формам собственных колебаний 
круглой пластины с закреплённым краем, аналогично работе [20]: 


и (г) = у и „Фо ии р)рар, (10) 
0 
где ф„ (г)=А, Л (К„г); К, — корни уравнения Л, (К„)=0,а А, =\№2/ 2 (К„) для тем. 
Учитывая линейность задачи, разложим функцию контактных напряжений (2) в ряд следу- 
ющего вида: 


Ч“) = Уи, а, (), 0</<1. (11) 


Здесь функции 4, (г) находятся из решения парного интегрального уравнения (9) с заданной пра- 


вой частью вида (10). 
Для построения решения рассмотрим следующее парное интегральное уравнение: 


Го (в) (в^) 2 (ега- 051+ 5, й бы, и 51, 
т (12) 
Го (в) 2 (аг)да=0, г>1. 


Из результатов работы [21] следует, что при построении решения задачи двухсторонне 
асимптотическим методом в качестве аппроксимации главной части трансформанты ядра можно 
использовать аналитическое выражение вида 


1, (и) = 4 ый (13) 


т тр 





где Ь,, и,, О — некоторые константы. 


В работе [22] получено аналитическое решение парного интегрального уравнения (12) в 
виде: 
26, (0)5 


Г : Вы (14) 


п 1-г Г 





Здесь введено обозначение: 





д 
Сумма в (14), содержащая функцию о ‚ соответствует решению уравнения (12) с пра- 
вой частью в виде ряда по функциям Бесселя. Функцию Ф„ (г) можно представить в следующем 
виде: 
ФЕ Л (г) р 


где и’ =А,; м’ =К». 


Таким образом, решение парного интегрального уравнения (9) с правой частью (10) имеет 
Вид: 
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а" (г) = 28175 [4 }(1- =)" Ус (мах | 
ат (г) =21*А,5 | [м (м„)Ф(г,Ж,) «У сич(там”) мы... 
1=0 
Здесь постоянные С;’ определяются из системы линейных алгебраических уравнений: 


>`с’а (ал „ВА“ +1 (0) 6 =0, К=1,2,...,М, 


№ 


У са АВВ КЕМ ТЕ 


1 
1= 


где 
а(а,6) = [азпа+ спа | (6? а ‚В(а,Ь) = [; (ла) а(а,Ь). 
Полученные значения контактных напряжений 4” (г) и заданную внешнюю нагрузку, при- 


ложенную к пластине, можно представить в виде следующих рядов: 


9» (г) = >. Ф, ( 7-91 р)рар, (15) 
ам (г), р, = [р()Ф, (р)рар. (16) 


Подставив разложения (10), (15), (16) в (8), придём к бесконечной системе линейных ал- 
гебраических уравнений для определения коэффициентов и/„ ‚ которую можно записать в следу- 
ющем виде: 


и” таки, Ел =р,К», т=0,1,2,...; а=-1, (17) 


=2п А, 5А, Ее (^к,)х(к,,к„) +5 С" х 2] 
п=0 


где ]=1,2,...; т=0,1,2,...; 





х(а,в) 5 [совать + (51а Ь)(а ь)" (2+5) +5)" 


Х (а,Ь) = (а° +52) (азпасозЬ + Бэт$ спа). 


В частности, 


= 21" /25А, Ё (0) кт мА, +51 Х (а, ^7,К .)| 


п= 
Е = тт] 1 (0) + са, ын(а,^') р 
м : 
 =21 124,5 | (Лк, Кл пк, + ло. ра. +м(а.^ 
Используя метод редукции, сведём решение бесконечной системы алгебраических уравне- 


ний (17) к решению системы линейных алгебраических уравнений следующего вида: 


М 
ИЕ, =р т = 0,1,2,...,М. (18) 


т! 
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После определения коэффициентов и/,, (т =0,1,...,М ) для фиксированного значения М 


из (18) находим контактные напряжения 


М 
9(г)= У и„9» (г), 0<г<1, (19) 
т=0 
и прогибы пластины 
М 
и (г) = У и„Ф„ (г), 0<г<1. (20) 
т=0 


Полученные формулы могут быть использованы для расчёта характеристик контактного 
взаимодействия многослойной пластины с основанием сложной структуры в зависимости от гра- 
ничных условий и характера нагрузки на пластину. Так как неоднородность пластины и многослой- 
ность в рамках рассмотренной модели влияет только на величину цилиндрической жёсткости 0, 
будем считать, что увеличение количества слоёв в пластине влияет на значение параметра Ди, 
соответственно, на значение параметра 5. 

Рассмотрим случай влияния на распределение контактных напряжений под пластиной рав- 
номерно распределённой нагрузки. В случае, когда края пластинки не закреплены по периметру, 
при вычислении контактных напряжений использовался результат работы [16]. Для закреплённых 
краёв использовалась формула (19). 

Пусть пластина лежит на мягком покрытии (В = 100), модуль Юнга которого меняется по 
закону: 


Е, (2) =Е,-Ф(2), Ф(2)=Ф, +(Ф -1).2. 
Показатель неоднородности ф, = сопз > 1 характеризует отношение модуля Юнга на поверхности 
покрытия (2 = 0) к модулю Юнга подложки. Ограничимся рассмотрением случая ф, = 3,5 , что со- 


ответствует сочетанию мягкого (алюминий, серебро, медь, свинец и т. д.) и жёсткого (железо, 
сталь, палладий, молибден и т. д.) металлов. 

При численной реализации было выявлено, что граничные условия для жёстких пластин 
практически не оказывают влияние на распределение контактных напряжений. Причём диапазоны 
параметров, для которых это справедливо, следующие: для Л = 4 — 5< 0,1, длял = 0,05 — 5< 0,07. 

На рис. 1, 2 построены графики контактных напряжений в случае гибких пластин для ука- 
занного выше закона изменения модуля Юнга с глубиной. Из графиков видно, что условия закреп- 
ления краёв пластины оказывают влияние на распределение контактных напряжений под пласти- 
ной. 

Заключение. Получено приближённое аналитическое решение задачи об изгибе круглой много- 
слойной пластины постоянной толщины, лежащей на упругом неоднородном основании сложной 


структуры. 


9(г) 










м $=4, 7120.05 


1.00- 
0.98- 
О Е =”) 
0.94- 


0.92 





0 0.10.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 г 











—— Закрепленные края —*— Свободные края 





Рис. 1. Распределение контактных напряжений под гибкой ($ = 4) пластиной 
при большой зоне контакта относительно толщины покрытия 
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Рис. 2. Распределение контактных напряжений под гибкой пластиной 
при малой зоне контакта относительно толщины покрытия 


Построенная математическая модель эффективна в широком диапазоне значений как гео- 


метрического (Л), так и физических (5, В) параметров задачи и может быть использована для рас- 
чётов неоднородных пластин Кирхгофа при различных граничных условиях на краю пластинки. 
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АХТ5УММЕТКТС ВЕМОТМС ОЕ СТВСУГАВ ЗАМОМЛСН РЕАТЕ ОМ ЕГА$ЗТТС РОУМОАТТОМ 
МЛТН СОМРЫСАТЕО $ТВОСТУВЕ* 


$. М. АмНомтСН, $. $. МоКом, А. У. МеКопуап 


Тре арргохтае с/озед-ютт зоиНоп {0 те ргоМет оп а сисшаг ти! ауег (5апаитсй) сопапЕ-/скпез$ рее 
Бепапд оп ап еазНс юипаайоп ий" а сотрсейеа $гисЁиге 5 оМатеад. ТпВе р/айе /5 БепЁ ипдег {Ве ахвуттес 
@УЕБИЕеА сад, апа ве юипаайоп геасйоп. Тйе е/а5Вс Гоипаавоп [5 а ип!оптиу йтедиаг п &/скпезс [ауег (соайпд) 
Базе оп а Вотодепеоиз Ва!-5расе (зибх{гае). Уоипд’5 тоди/из уаие аЁ те 1тегасе о! {те соаНпд апа те 
5ирзтее Ваз а чдтЙсапЕ [еар. Тио аТегепЕ сазез о! пе Боипаагу сопюопз аге сопз!аегеа Гог те ре: йхеа 
апа ее еаде сопопв. Тре сопзгисеа арргохмтае апа!уйса! зоиНоп {0 те рго ет /5 еТесьуе Гог а илае гапде 
ОЕБоЁ деотеё с (те пвотодепеои$ [ауег Ёи/скпе55 апа {пе р/аЁе гаи) апа рпуз/са/ рагатеег (р/айе ПехИКу, 
апа е/аНс ргорегие$ оГ те соайпд апа {Те зибзгае). Тпе сотасЕ ргоШет /5 гедисеа 0 {пе зучет оЕте педго- 
ЯЙТегепйа/ едиайоп зо/иНоп {&гоидв {пе птедга! гапзюЮгтайоп тетпоа. Тйе оБбатеа юпти/а$ сап Бе изей Юг 
сасшайта те сотщасЕ и®егасйоп спагасетсНс$ Бериееп а ти вуег ре апа а Гоипаайоп ив а сотр/ех $#ис- 
{ге тп уапоиз сазез о {те Боипдагу сопопз, апа иапоис [оа45 арр/еа ю #е рае. 

Кеуигогаб5: ппотодепеоиз таейа(, запаиисй р/ае, Гипсвопа/у-дгадеа соайпа, ахзуттеё ргоМет, апа/уйса! 
теосб, арргохтае апа/У#са! 5о/ийоп. 





* Тне гезеагсй 15 зиррокед Бу ВЕЕТ дгапе по. 13-08-90916-то! {п_пг. 
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